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On the sizes of z > 2 Damped Lya Absorbing Galaxies* 
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ABSTRACT 

Recently, the number of detected galaxy counterparts of z > 2 Damped Lyman-a Ab- 
sorbers in QSO spectra has increased substantially so that we today have a sample of 10 de- 
tections. M0ller et al. in 2004 made the prediction, based on a hint of a luminosity-metallicity 
relation for DLAs, that Hi size should increase with increasing metallicity. In this paper we 
investigate the distribution of impact parameter and metallicity that would result from the 
correlation between galaxy size and metallicity. We compare our observations with simu- 
lated data sets given the relation of size and metallicity. The observed sample presented here 
supports the metallicity-size prediction: The present sample of DLA galaxies is consistent 
with the model distribution. Our data also show a strong relation between impact parameter 
and column density of Hi. We furthermore compare the observations with several numerical 
simulations and demonstrate that the observations support a scenario where the relation be- 
tween size and metallicity is driven by feedback mechanisms controlling the star-formation 
efficiency and outflow of enriched gas. 

Key words: galaxies: formation - galaxies: high-redshift - galaxies: ISM - quasars: absorp- 
tion lines - quasars: individual: SDSS J 091826. 16+163609.0, SDSS J 033854.77-000521.0, 
PKS0458— 02 - cosmology: observations 



1 INTRODUCTION 

The size of high redshift galaxies, its correlation with other proper- 
ties like mass or luminosity, and its e volution with redshift i s a fun- 
damental quantity to measure (e.g., iRees & OstrikeJ[l977h . QSO 
absorption-line studies provide an interesting way to address these 
relations as the method is s ensitive to galaxies ove r a large range 
of masses and s izes (e.g., Nagamin e et al.l [20071 : IPontzen et all 
2OO8: ICen|[2O10h and because it is complementary to the infor- 
mation gathered from studying the stellar light of galaxies (e.g., 
iFerguson et ai]|2004l) . 

Damped Lyman-a Absorbers (DLAs) are believed to be the 
absorption s ignatu re of neutral interstellar gas in galaxies (see 
Wolfe et al] J2005h for a review). When the first sample of DLAs 
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was collected it was found that the total cross-section for DLA ab- 
sorption was about five times higher than what would be expected 
from the density and sizes of local spiral galaxi es ( Wolfe et al.l 
198d see also Mol ler & Warren 1998). IWolfe et al.ldl986l) and also 



Smith et all J1986I) put forward the view that DLAs could be form- 
ing disks with radii ~ v5 times larger than local spirals. In this 
view galaxies were larger in the past, which would be at odds with 
the pr evailing hiera rchical model of galaxy formation. On the other 
hand, iTvsonl J 1988h argued that DLAs could be caused by a pop- 
ulation of many small (much smaller than local spirals), gas-ric h 
dwarf galaxies (see also lHaehnelt et all 1 9981 : [Ledoux et al.l ll998). 

Detecting the absorbing galaxies in emission is the most ob- 
vious way to determine the sizes of DLAs from the distribu- 
tion of impact parameters relative to the background QSOs. An- 
other method is to measure the extent of 21 cm or X-ray ab- 
sorption against ext ended background sources (iBriggs etalJI 19891 : 
Diik straetai1l2005h . but such extended background sources are 
currently rare. Also, informatio n on sizes can be i nferred from 
DLA s towards lensed QSOs dSmette et all 1 19951 : lEllison et af] 
|2004 120071: ICooke et alj fcoifA but such systems are also rare. 
Given the proximity of the bright QSO to the line of sight to the 
absorbing galaxy the first approach adopted was to search for Lya 
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emission from the absorbing galaxy as t he DLA itself he re re- 
moves the light from the background QSO (Foltz et al. 1986). Two 
observing strategies were adopted reflecting two different ideas 
about the nature of DLAs. Under the assumption that the Hi-extent 
of DLA galaxies is small their impact parameters relative to the 
QSO will also be small. Therefore a long-slit centered on the QSO 
should have a high probability of also covering the DLA (e.g. 
iHunstead et alj [l990. and references therein). If on the other hand, 
DLAs arise from large Hi disks they may have large impact param- 
eters relative to the QSO and narrow-band imaging will be a better 
approach as it eliminates the risk of missi ng the object on the slit 
(Smith et al. 1986; M0 ller & Warrenl 19931) . Later , spectroscopy us- 
ing integral field units has also been attemp ted (Chris tensen et al.l 
l2007l|Peroux et al.ll201ll ; lBouche et ai]|2012l) . 

In the period 1986-2010 only two galaxy counterparts of 
bona-fide intervening DLAs at z > 2 (log A<Hi/cm~ 2 > 20.3, 
Zabs « z em an d AV > 6000km s _1 ) were detected (see 
M0lle r et alj |2004|) despite many more systems were observed. 
This meager detection rate has gradually been understood to re- 
flect the fact that DLA galaxies, due to their cross-section se- 
lection, trace the faint end of the galax y luminosity func tion 
(Fynbo et al. 2008, and references therein). lM0ller et al.l(l2004h and 
iLedoux et al] j2006l) found tentative evidence that DLA galaxies 
obey a metallicity-luminosity relation similar to local galaxies and 
suggested that a survey specifically targeting high metallicity DLAs 
should lead to a higher success-rate in searches for galaxy counter- 
parts. Indeed, this strategy has more than doubled the number of 
detections of z > 2 DLA galaxy counterp arts in the literature in 
the last two years dFvnbo et alj201Ctl201 II) . 

The objective of this work is to test the correlation be- 
tween metall i city a nd physical size of DLA galaxies predicted by 
M0lle r~et alj d2004l) . Thanks to the recent increase in number of 
detections of emission from DLAs we can now analyze how the 
distribution of impact parameters and their corresponding metal- 
lic ities match up wit h the expected distribution from the model 
of iFvnbo et al. liooi), based on the assumption of the previously 
mentioned correlation. 



2 OBSERVATIONS AND SAMPLE SELECTION 

Q 0338-0005 was observed on November 9 2010 with the X- 
shooter spectrograph mounted on UT2 of ESOs Very Large Tele- 
scope follo wing the same strat egy as that described in lFvnbo et al.1 
(2010) and Fynb o et alj J201 l|), i.e. securing 1 hr of integration at 
three position angles 0° and ±60°. For details on t he reduction 
we ref er to the description in IFvnbo et"aL I fcOlOb and lFvnbo etall 
fcOllh . 

We detect the galaxy counterpart of the DLA in the trough of 
the damped Lya line in all three spectra (PAs of 0°, ±60°). The 
lower panel of fig. Q] shows the stacked 2D spectrum around Lya 
where the Lya emission line can be seen in th e bottom of the DLA 
tro ugh. Using the tria ngulation described in M0 ller et al. 1 d2004l) 
and lFvnbo etall J20 id) we infer an impact parameter of 0.49 ±0. 12 
arcsec. The metallicity is determined to be [M/H] = —1.25 ± 0.10 
based on low-ionization absorption lines from Sill in a high res- 
olution UVES spectrum of the QSO (Ledoux, private commu- 
nication) and on the damped Lya line from which we measure 
log A^m/cm" 2 = 21.05 ± 0.05. 

The target PKS0458-02 was observed on November 24 2009 
with X-shooter. The slit was placed at a position angle of —60.4° 
EofN with the purpose of measuring the precise impact parameter 
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Figure 1. Two dimensional spectra from X-shooter of PKS0458— 02 (top) 
and Q0338— 0005 (bottom) shifted to rest-frame wavelengths. In the centre 
of each spectrum the DLA troughs are seen with Lya emission from the ab- 
sorbing galaxies at spatial offsets of 0.31±0.04 arcsec and 0.49±0.12 arc- 
sec, respectively. The spectra has been smoothed for illustrative purposes. 



of the z = 2.04 DLA galaxy. In the X-shooter spectrum (shown 
in the top panel of fig. [T) we measure an impact p arameter of 
0.31± 0.04 arcsec. We have combined the data from M0ll er et al.l 
d2004l) with our new measurement of the impact parameter. 

The galaxy counterp art of the z = 2.58 DLA towards 
Q0918+1636, reported in IFvnbo etail d201ll) . has here been ob- 
served in a ground based imaging study using Point Spread Func- 
tion (PSF) subtraction to secure a more precise measure of the im- 
pact parameter. The field was observed with the Nordic Optical 
Telescope (NOT) using the ALFOSC instrument in service mode. 
The target was observed in the R-band over the nights March 27, 
29 and April 9 201 1, giving a total of 16 380 sec exposure. The im- 
ages were reduced and combined using standard methods for CCD 
imaging data. 

The field around Q0918+1636 is shown in the left panel of 
Fig. [2] The figure shows a 20x20 arcsec 2 region from the com- 
bined NOT image. In the right panel we show the field after PSF- 
subtraction and smoothing with a 3x3 pixel box-car filter. Two 
sources, both with i?-band magnitudes of approximately 25, are 
detected at impact parameters of 2.0 arcsec and 3.5 arcsec. Closer 
to the position of the QSO there are residuals from the PSF subtrac- 
tion. With dashed line s we show the pos itions of the three X-shooter 
slit positions used bv lFvnbo et al. lionl). As seen, the continuum 
source at 2.0 arcsec seen in Fig. |2]coincides both in position angle 
and impact parameter with the emission line source detected in the 
X-shooter spectrum and must hence be the continuum counterpart 
of the z — 2.58 DLA galaxy. The source at 3.5 arcsec falls outside 
of the X-shooter slits and it is on the basis of the data in hand not 
possible to establish if this source may be related to t he z = 2.41 
DLA also detected in the spectrum of Q0918+1636 dFvnbo et"al] 
l201ll) . 

2.1 Sample selection 

In this paper we base our analysis on data combined from the liter- 
ature along with our new data for three systems of DLA galaxies; 
Q0338-0035, Q0918+1636 and PKS0458-02. In addition to the 
systems in the literature we include a con firmed candidate from an 
expanded sample based on the methods in lChristensen et al. I d2007l) 
that has been confirmed by long-slit spectroscopy (Christensen et 
al., in preparation). We note, however, that no other candidates 
from that work has been confirmed. Bunker et al. in preparation 
also have a detection o f Lya emission from a DLA galaxy using 
long-slit spectroscopy (Weatherlev et al. 2005, their table 1). The 
redshift, metallicity and impact parameter of the systems used in 
our analysis are shown in Table [TJ The z a bs ~ z em DLA towards 
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Figure 2. R-band image of the 20 X 20 arcsec 2 field of Q0918+1636. 
The left panel shows the image before PSF-subtraction and the right panel 
after PSF-subtraction and smoothing with a 3 X 3 pixel box-car filte r. Also 
plotted are the three X-shooter slits used in the iFvnbo et alj (2011) study. 
The galaxy counterpart of the z = 2.58 DLA detected in the X-shooter 
spectrum taken at PA = +60° is indicated with an arrow. Also seen is 
another galaxy at a PA of about 163° and an impact parameter of 3.5 arcsec. 



Table 1. z > 2 DLAs with identified galaxy counterparts used in this study. 



References: (1) Irvfeller et alj J2004I); 
IFvnbo et allfeOllI): (4)lBouche et al j 20K 
et al. i n prep.; (7)lNoterdaeme*et*aU i 20n 



(2) IFvnbo et alj feOld) ; (3) 
(5) this work; (6) Christensen 
(8) L edoux, priv. comm.; (9 ) 



IWolfeetalJb008l) ; (10) Bunker, priv. comm.; (11) Prochask a'etai] J2003l) . 



PKS0528— 250 has been shown to be unrelated t o the nearby QSO 
and hence likely is similar to intervening DLAs (M0ll er & Warrenl 
ll993l ; lM0ller & Warrenll998h . 



3 THE METALLICITY-SIZE CORRELATION 

The re lation between metallicity and size modeled by F vnbo et al.1 
is based on the assumption that DLAs seen in the spectra 
of background QSOs are caused by galaxies comparable to low lu- 
minosity Lyman Break Galaxies. In this model the galaxies are ap- 
proximated by randomly inclined disks, each one assigned a size, 
metallicity and metallicity gradient based on its luminosity. The 
different relations between the quantities listed above are all based 
on local observationally motivated correlations for which indica- 
tions of their validity at higher redshift have been seen. The more 
luminous galaxies have bigger disks, higher metallicities and shal- 
lower metallicity gradients. The luminosity L is then drawn from 
the LBG luminosity function weighted by Hi cross section. An im- 
pact parameter is also assigned again with the luminosity depen- 
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Figure 4. Distribution of (top) slope vs. intercept and (bottom) slope vs. 
scatter from the 4000 simulated data sets (grey points). The solid and 
dashed lines represent the 1 a and 2<r confidence contours, respectively. The 
black cross shows the result from our fit to the data. 



dent cross section of Hi as weighting factor. The random sight line 
is then given a metallicity from the before mentioned relations con- 
necting [M/H] to L and the radial d ependence of metallicity. For 
further details see lFvnbo et al. l (E008h . 



3.1 Testing the model 

Though the sample we present here is currently the largest at hand, 
we still only have ten points. Therefore we need a simple statistical 
estimator to describe the data. We have chosen a two-parameter 
fit of the form: log(b) = a + /3-fM/H]. Since the internal scatter 
in the data is larger than the individual measurement errors we fit 
a straight line to the measured data by least-squares minimization 
without weights. 

We then use the model described above to generate simulated 
data sets to which we can compare the measured data. Because the 
measured data points are selected in different ways we cannot sim- 
ply draw ten random points from the model. In order to make a 
valid comparison, we choose model points around the same metal- 
licities as the measured points. 

The simulated data points have metallicity taken randomly 
from a gaussian distribution around each measured point using the 
measurement error as the width of the distribution. We then take 
the distribution of impact parameters from the model at the given 
metallicity and randomly pick an impact parameter. We do this for 
all ten points to simulate a complete data set 4 000 times. Every set 
of points is then fitted by a straight line to give the intercept, slope 
and scatter. 

The 4 000 simulated data sets are shown in the right panel of 
fig. [3] as blue points along with the measured data points in black. 
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Figure 3. Impact parameters plotted against Hi column density (left) and metallicity (right) for z > 2 DLAs with securely identified galaxy counterparts, see 
table [T] The observed impact parameters have been converted to kpc at redshift 2 = 3 instead of arcsec for easier comparison. Th e blue points in th e right 
panel show the simulated distribution of impact param eters as a function of metallicity for DLA galaxies at z = 3 from the model in Fynbo et al. (2008). The 
red points in the left panel show model points from lPontzen e t al. 1 2008). 



To make the plot less crowded only every fifth simulated data set 
has been shown. The left panel of fig.[3]shows impact parameter vs. 
neutral hydrogen column density for our data in big black points. 
The re d small points in the panel are results from IPontzen et al.l 
d2008h . 

In fig. [4] we show the results of the fits to the data, both 
observed and simulated. The figure shows the la (solid) and 
2a (dashed) contours of the distribution (grey points) of fitted 
slope and intercept and slope and scatter. The black cross in each 
panel point represents the best fitting values to the data: a = 
1.15 ± 0.70 log(kpc), P = 0.41 ± 0.67 log(kpc) dex -1 with 
a scatter around the fit of 0.26. The mean values of the intercept 
and slope from the simulated data are 0.90 ± 0.28 log(kpc) and 
0.41 ± 0.26 log(kpc) dex -1 , respectively, with a mean scatter of 
0.27. 



4 DISCUSSION 

The sample of spectroscopically confirmed z > 2 DLA galaxies 
is now large enough to start determining statistical properties of 
DLA galaxies. Here we have presented a new detection, namely of 
the DLA galaxy counterpart of the relatively metal poor DLA to- 
wards Q 0338-0005 and have detected the z = 2.58 DLA galaxy 
towards Q 0918+1636 in the continuum. Based on these and previ- 
ous detections we have analyzed the distribution of impact parame- 
ters vs. metallicity to study the underlying correlation between size 
and metallicity. The results presented in fig. |4] show that our data 
are consistent with the model at the 1.2<j level for slope vs. inter- 
cept. Our data is well within the lcr-contour in slope vs. scatter. The 
agreement is also seen in the distribution shown in the right panel of 
fig. [3] There is, however, a single data point which is significantly 
outside the model distribution (see fig.|3] right panel, top part). This 



can be explained as a result of the simplicity of the model. The cor- 
relations that go into the model all have internal scatter. Including 
this scatter will result in a larger spread in physical size, hence mak- 
ing the upper edge of the model distribution less tight, and thereby 
including the outlying point. 

An important point to note is that four of the data points 
have specifi cally been chosen to have high metallicity ( Fynbo et al. 
|2010L 1201 lh . When comparing the observations with the model it 
is therefore important to remember that the observed sample is bi- 
ased towards high-metallicity DLA galaxies that in our model are 
expected at the largest impact parameters. At the same time we 
may be l oosing points at the high-metallic ity end because of dust 
bias (e.g. JPei et all 1 9991 ; iKhareet alj20llh and due to the fact that 
our method using th ree long-slits is les s sensitive to high impact 
parameter DLAs, see lFvnbo et al. I d2010l) for a detailed description 
of the method. Nevertheless, the data on high-luminosity DLAs are 
consistent with the picture in which gas at large impact parameters 
of tens of kpc cause DLA absorption. These galaxies seem to fol- 
low relations between size, luminosity and metallicity comparable 
to the relations we see for disk galaxies in the local universe. 

From hydrodynamic simulations it has also be en investi- 
gated in which kind of systems DLAs originate (e.g., Cenl 20 lOl : 



IPontzen et al.l2008l;lRazoumov et alj200d : lGardner et al.l20oTh . In 
the simulation by ICenl (2010|) DLAs at z = 3 — 4 are dominated 
by low metallicity systems at large impact parameters, typically 
b — 20 — 30 kpc. In this simulation the contribution from gaseous 
disks to the DLA incidences is very small. Instead most DLAs arise 
in filament ary gas. The impact parameters at low metallicity from 
ICenl ( 1201 Of) are much larger than what we observe. The only metal 
poor DLA is found at a n impact pa r ameter of 2.5 kpc, henc e not 
supporting the model bv lCenl bOlOh . lRazoumov et al.l I2OO6I) also 
find a large number of DLAs at high impact parameter from fil- 
amentary gas around large structures. But they find that the im- 
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portance of these filamentary DLAs decreases as they increase the 
reso lution, 

IPontzen et alj J2008h are able to reproduce the properties of 
DLAs well with their simulation, which contains very detailed 
feedback mechanisms. From their simulation they see positive cor- 
relations between the DLA cross-section and the mass of the halo 
containing the DLA and between the metallicity and the halo mass 
(see their fig. 4 and fig. 18). From the correlation between ctdla 
and M v i r it is possible to infer a limit on the impact parameters 
of the corresponding DLAs by assuming a spherical gas disk with 
a radius given by: J?dla oc yf ctdla- The largest possible cross- 
section found in the most massive haloes is <r ma x ~ 10 3 kpc 2 . 
Hence, the impact parameters for DLAs around z ~ 3 in their sim- 
ulation are expected to be b < 30kpc. The correlation of [M/H] 
and halo mass suggests that the metal poor systems are to be found 
in smaller haloes and therefore at smaller impact parameters. The 
opposite is true for metal rich DLAs. This agrees well with the data 
presented here in Fig. [3] 

We see a fairly strong anti-correlation between log(iVHi) and 
impact parameter (Spearman rank = —0 .6), see left panel i n fig.[3] 
A similar trend is seen in studies by llVfoller & Warrenl dl998l) ; 
Moni eret alj d2009l) and bv lPeroux et al.M201l|). however, we find 
a tighter and steeper correlation in our data. iGardner et al.l |200 1 ) 
also find a clear anti-correlation between impact parameter and 
log(NHi) in their simulation of DLAs in different cosmologies, 
including ACDM. They too find a shallower relation for DLAs 
at z = 3, This is most probably a result of the very simpli- 
fied algorithms to ha ndle cooling and feedback in their simulation. 
IPontzen et al.1 {2008) find a clear relation between log(NHi) and b 
(red points of fig. [3} that coincides very well with our observed 
distribution as seen in the left pane l of fig. [3] The good agreement 
between data and the simulation of Pontze n et al.1 d2008h indicates 
that the feedback mechanisms in their simulation that control star 
formation and outflow of enriched gas are responsible for the cor- 
relation between size and metallicity. 

In conclusion, our analysis shows that our model well predicts 
the statistical properties of the observed sample without adjustment 
of any parameters, and as such provides strong support for the un- 
derlying size-metallicity relation. The nature of the metal poor pop- 
ulation of DLA galaxies is, however, still unknown, as these DLA 
counterparts are much more difficult to observe. 

Furthermore, this sample of spectroscopically confirmed DLA 
galaxies serves as constraints for future simulations that include 
modeling of DLAs. 
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